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Abstract 

In this work, we have studied theoretically the effects of gold adsorption on the Al(OOl) surface, 
using ah initio pseudo-potential method in the framework of the density functional theory. Having 
found the hollow sites at the Al(OOl) surface as the most preferred adsorption sites, we have 
investigated the effects of the Au adsorption with different coverages (0=0.11, 0.25, 0.50, 0.75, 
1.00 ML) on the geometry, adsorption energy, surface dipole moment, and the work-function of 
the Al(001) surface. The results show that, even though the work-function of the Al substrate 
increases with the Au coverage, the surface dipole moment decreases with the changes in coverage 
from = 0.11 ML to = 0.25 ML. We have explained this behavior by analyzing the electronic and 
ionic charge distributions. Furthermore, by studying the diffusion of Au atoms in to the substrate, 
we have shown that at room temperature the diffusion rate of Au atoms in to the substrate is 
negligible but, increasing the temperature to about 200° C the Au atoms significantly diffuse in to 
the substrate, in agreement with the experiment. 

PACS numbers: 68.43.Bc, 81.65.Mq, 68.47.Gh 
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I. INTRODUCTION 



Modification of the Al surface properties by deposition of different kinds of atoms has been 
studied both experimentallyi^ 1 ^^ and theoretically.-*^ 1 ^ Specifically, the deposition of 
Au atoms on the Al(OOl) surface, which is the subject of our present work, has been studied 
experimentally^ and to our knowledge, there is no report in the literature on theoretical 
work. The analysis of the experiment conducted by Egelhoff 6 had shown that for sufficiently 
high coverages (around 10 ML), the Au atoms had been deposited on the Al(001) surface 
in the form of clusters which is different from the layer-by-layer deposition. 

In this work, we have studied the effects of deposition of gold atoms with low-coverages 
(0 = 0.11, 0.25, 0.50, 0.75, 1.00 ML) on the Al(001) surface properties. At these low 
coverages, the clustering effect would be negligible. For this purpose, we have carried out 
first-principles calculations within the density functional theory (DFT).— ^ Our calculations 
for the clean Al(001) surface shows outward relaxations for the interlayer spacings, consis- 
tent with the other ab initio calculations^^^^ 1 ^^^ and the latest experiment .-^ In the 
Au/Al(001) system, we have shown that the work-function, average binding energy, geome- 
tries, and induced surface dipole moments of the Al(001) surface changes with increasing the 
surface density of Au adatoms. Since the electronegativity of Au atom is higher than that of 



Al atom (2.4 for Au and 1.5 for Al, See Ref. 



22l ). adding Au atoms on the Al surface, pulls 



more electrons from the substrate and hence increases the magnitude of the surface dipole 
moment which, in turn, increases the work-function. However, by increasing the surface 
density of Au adatoms, the induced dipole moment shows a decrease from = 0.11 ML 
to = 0.25 ML with a sharp dip at 0=0.25 ML, a maximum at 0=0.50 ML, and a slow 
decrease up to 0=1.00 ML. We have explained this behavior by analyzing the distribution 
of space charge of electrons and ions in the slab. The organization of this paper is as follows. 
In Sec. Hi] we explain the detailed calculational methods. Sec. [Till is devoted to the results 
and discussions and finally, we conclude this work in Sec. HVl 

II. CALCULATION DETAILS 

We have used the Quantum-ESPRESSO code^ for the electronic structure calculations 
of the bulk Al, Al(001) surface, and Au/Al(001) system. The calculations are based on the 
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DFT and the self-consistent solution of the Kohn-Sham (KS)^ equations. For the exchange- 
correlation, we have used the GGA of Perdew et aM^\ and for the valence electrons of the Al 
and Au atoms, we have used the scalar-relativistic ultra-soft 2 ^ pseudo-potentials.™ For the 
Al pseudo-potential, the electronic configuration 3s 2 3p l with the same cut-off, r c = 1.76 a.u., 
for both orbitals; and for the Au pseudo-potential, the electronic configuration 6s 1 6p°' 5 5d 9 ' 5 
with respective cut-offs r c = 2.4, 2.6, 2.0 a.u. were used. 

To calculate the bulk properties of fee Al, we have tested different k-point samplings and 
energy cut-offs for the plane wave basis, and have performed the Brillouin-zone integrations 
using the Methfessel-Paxton smearing 27 with 0.05 Ry of broadening for appropriate set of 
(20 x 20 x 20) Monkhorst-Pack gridM and a 38 Ry cut-off for the plane wave basis set. 
Using these parameters, we have solved the self-consistent KS equations for different lattice 
constant values and obtained the equilibrium lattice constants as those minimizing the total 
energies. We have considered the atoms in their equilibrium state, when the forces on the 
atoms were converged to within 1 mRy/a.u. 

The Al(001) surface is modeled by a supercell composed of a 7 atomic layers with a 
(lxl) surface unit cell (with the surface-cell size taken from the bulk calculation) and a 
vacuum region of 15 empty layers. The kinetic energy cut-off for the plane wave basis and 
the k-point mesh were chosen as 38 Ry and (20 x 20 x 1), respectively. 

For the Au/Al(001) system, using Eq. (JT]), we have first found the hollow sites as the 
preferred adsorption sites, and then a 7-layered Al slab was taken as substrate. The Au 
atoms were positioned above the preferred adsorption sites on the Al(001) surface, and a 
vacuum region of 13 empty layers was taken. The surface structures of (3 x 3)-Au, (2 x 2)-Au, 
(2 x 1)-Au, (2 x 2)-3Au, and (1 x 1)-Au were chosen for coverages 6=0.11, 0.25, 0.50, 0.75, 
and 1.00 ML, respectively. In order to preserve the inversion symmetry, we have positioned 
the Au atoms on both sides of the slab. For Au adsorption on the Al(001) surface, the 
average binding energy per Au atom is defined as 

£Au/Al(001)^ 3_[£Au/A1(001)(q) 

Navl 

-(N An E Au + mi m 2 E A1 ^)}, 

in which N Au is the total number of symmetrically adsorbed Au atoms on the two sides of 
the slab (N Au /2 on each side) and E An is the total energy of an isolated Au atom. _E A1 (° 01 ) is 
the total energy of a 7-layered Al supercell having a (lxl) surface unit cell with one Al atom 

3 



at each layer, and #a u /ai(ooi) (q) j s the total energy of a supercell composed of a 7-layered Al 
substrate and Au adatoms with coverage B and the surface structure of Al(001)-(mi x WI2)- 
(iA^Au)Au. For example, in 6 = 0.75 ML, £ a u/ai(ooi)(@) is the total energy Q f a sup ercell 
composed of a 7-layered Al substrate with 4 Al atoms on each layer and 3 Au atoms at 
each side of the slab. The surface structure in this case is Al(001)-(2 x 2)-3Au. The binding 
energy, Eq. (0Q), is the energy that is gained upon adsorption of a free Au atom on the Al(001) 
surface. In this definition, a positive number indicates that the adsorption is exothermic 
with respect to the free Au atom; and a negative number indicates that it is endothermic 
(i.e., adsorption does not take place). To evaluate the energy of an isolated Au atom, E\ u , 
we have solved the self-consistent KS equations for an Au atom in a cubic supercell of 20 a.u. 
side, a (7 x 7 x 7) k-point mesh for the Brillouin-zone integrations, and 38 Ry for the cut-off 
energy of the plane- wave basis. This choice of k-point mesh was adopted so that the k-point 
density becomes the same as in the Au bulk calculations. However, the result obtained using 
the Gamma point alone was the same to within 2 decimals. On the other hand, choosing 7- 
and 9-layered Al slabs and adsorption of 1 ML of Au atoms in each case, the difference in 
the binding energies was within 10 -3 eV, and the difference in the equilibrium distance of 
Au atom from the Al surface was less than 0.01 A. With these testings, we have chosen a 
7-layered Al slab as the substrate in our calculations; a (20 x 20 x 1) k-mesh, and 38 Ry cut- 
off energy for studying the Au adsorption with (lxl) surface structure which corresponds 
to G=l ML. In order to preserve the same k-point density for different coverages, we have 
used (7x7x1) mesh for (3 x 3)-Au (6 = 0.11 ML), (10 x 10 x 1) mesh for (2 x 2)-Au 
(0 = 0.25 ML), (10 x 20 x 1) mesh for (2 x 1)-Au (6 = 0.50 ML), and (10 x 10 x 1) mesh 
for (2 x 2)-3Au (0 = 0.75 ML) surface structures. 

In our supercell calculations, to calculate the work-function of a given surface, which 
is defined as the minimum energy needed to remove an electron from a crystal to a point 
outside that surface, we find the difference between the electrostatic potential in the middle 
of the vacuum region and the Fermi energy of the slab. For more accurate results, we have 
considered sufficiently large vacuum region to prevent the overlap of the wave functions of 
the neighboring slabs. 

In order to analyze the bonding behavior for different coverages, we have considered the 
difference electron density 2 ^, An(r), which is sensitive to the changes in the charge density, 
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and is denned as 

An(r) = n Au / A1 (° 01 )(r) - n A1 (° 01 )(r) - n Au (r), (2) 

in which ri AU//A1 ^ 001 ^ (r) is the total electron density of the substrate-adsorbate system, 
n A1 ( 001 )(r) is the total electron density of clean Al(OOl) slab, and n Au (r) is the total electron 
density of the isolated Au adlayer. 

To evaluate the induced surface dipole moment due to Au adsorption, we adopt two 
different methods. In method 1, we use the Helmholtz equation (in Debye units),— which 
is defined as 

AAW 

in which A is the area in A 2 per (lxl) surface unit cell, AW is the work-function change 
in electron volts, and is the coverage. In method 2, we resort to the self-consistent 
charge-redistributions and directly use the integration 

/■Zvac 

/x ind = Ap z = dz z [p 2 (z) - pi(z)} , (4) 

Jz=0 

in which p 2 {z) is the integral of the self-consistent adsorbate-substrate charge density over 
the surface unit cell (m x n) divided by mn; and pi{z) is the same average, but here of 
the superposition of the self-consistent charge density of the isolated substrate and the self- 
consistent charge density of the isolated Au adlayer. The limits of the integral, z = and 
z VSLC are, respectively, the position of the central layer of the substrate and the distance of 
the middle of the vacuum region from the central layer in the adsorbate-substrate system. 
Here, the position of the central layer of the isolated substrate is identified with the position 
of the central layer of the adsorbate-substrate system and the position of the isolated Au 
adlayer is identified with the position of the Au adlayer in the adsorbate-substrate system. 
Therefore, because of fixed positions, the contributions of the ion charges of the Au adlayer 
and the central layer (in fact, in our case, the central as well as its first neighbor layers) in 
Eq. (J3J) cancel out. However, there is no such cancelation in the electronic charge density. 
In other words, in our 7-layered substrate, only the topmost two layers of the substrate 
contribute in the ionic redistribution. To apply Eq. (j4]), we use the discretized form 

» i 

in which Zi is the distance of the slice of thickness 5z = 0.01 A from the slab center and Aqi 
is the change in the slice electronic charge. The z'a and Zj are the ionic positions after and 



TABLE I: Lattice constant, bulk modulus, and cohesive energy of Al. 





present work 


LAPW a 


LAPW 6 


LAPW C 


PP-PW d 


PP-PW 6 


Experiment 


a (A) 


4.04 


4.04 


4.10 


4.04 


4.06 


4.05 


4.05±0.10 / 


B (GPa) 


73.8 


75 


73 


78 


75 


79 


72.2 s 


£coh (eV/atom) 


3.48 


3.65 




3.60 


3.52 


3.52 


3.399 
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a Ref. m 
b Ref. E- i 
c Ref. 
d Ref. 
e Ref. 

^Ref. 
s Ref. 137 
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before the Au adsorption, respectively. Here, we have taken the value of q = +3e for the Al 
ionic charge. 

III. RESULTS AND DISCUSSION 

In order to study the adsorption of Au atoms on Al(001) surface, we had to consider first 
the bulk Al, clean Al(001), and the isolated Au atom. 



A. Bulk Al and clean Al(OOl) surface 

Our calculated lattice constant and bulk modulus for the bulk Al is compared with 
experimental value and other theoretical results in Table. [H As is seen, our results are in 
good agreement with other theoretical works and experiment. 

In Al(001) surface calculation, we have obtained the relaxed positions of the first and 
second outermost layers, and using the value of the bulk spacing, we have calculated the 
interlayer relaxations, Ay = (d^ — d )/do x 100%, between layers i and j with respect to 
the bulk spacing (g?o=2.02 A). Our results show that both interlayer spacings, di 2 and g?23, 
have outward relaxations consistent with the latest LEED experimental result and other ab 
initio calculations. The values are presented in Table [Til Although our relaxation results for 
7-layered slab are smaller than the experimental and the 23-layered LAPW results (Ref. l2ll). 
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TABLE II: Relaxations of clean Al(OOl) surface. For comparison, other ah initio results (only 
GGA) as well as the experimental result are included. NL indicates the number N of atomic 
layers used in the slab, pp means pseudo-potential method. 





7 L-pp a 


7 L-pp 6 


7 L-LAPW c 


23 L-LAPW d 


Experiment e 


A12 


+0.91 


+0.5 


+1.477 


+1.51 


+2.0+0.8 


A23 


+0.26 


-0.3 


-0.057 


+0.42 


+1.2+0.7 



"Present, work 
6 Ref. M 



c Ref. 20. Calculations up to 17 L has been also considered for which A12 = +1.598 and A23 = +0.436 
^Ref. 
e Ref. 



but, in contrast to the results of Ref. 
the experiment. 



19| and Ref. |20j, the trend is in good agreement with 



B. Proper adsorption site on Al(OOl) surface 

In order to find the proper adsorption site of an Au atom on the Al(001) surface, we have 
considered a (2 x 2) surface unit cell for the coverage G=0.25 ML, and using Eq. (TO , we have 
obtained the adsorption energies of +2.58 eV/atom, +3.09 eV/atom, and +3.73 eV/atom 
for the top, bridge, and hollow sites, respectively. Since all these values are positive, there 
would be gain of energy in the deposition at each of the three sites. However, since the 
resulting values are quite distinct, we choose, confidently, the hollow sites as the preferred 
adsorption sites without any need to repeat the process for larger surface unit cells. We 
therefore, consider the adsorption at hollow sites for all coverages. In Fig. [H we have shown 
the side and top views of the surface structures (3x3)-Au, (2x2)-Au, (2xl)-Au, (2x2)-3Au, 
and (lxl)-Au, which correspond to the coverages of B=0.11, 0.25, 0.50, 0.75, and 1.00 ML, 
respectively. The black balls and grey balls represent Au and Al atoms, respectively. 

To study the diffusion of Au atoms into the Al substrate, we have considered a super cell 
with surface structure of (2 x 2) and put an Au atom at the hollow site under the first 
layer. Then we allowed the atoms relax until the force on each atom becomes less than 1 
mRy/a.u. The result showed that the gold atom has pushed the Al atom on its top toward 
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Coverage 0.11 ML 0.25 ML 0.5 ML 0.75 ML 1 ML 

(3x3)-Au (2x2)-Au (2xl)-Au (2x2)-3Au (lxl)-Au 



FIG. 1: The side and top views of the supercells used in the study of Au atoms on the Al(OOl) 
surface. The black and grey balls represent the Au and Al atoms, respectively. The supercells 
from left to right correspond to the coverages 0=0.11, 0.25, 0.50, 0.75, and 1.00 ML, respectively. 
Only the half of the Al substrate including the central layers are shown. The Al atoms labeled by 
"C" and "D", and the Au atoms labeled by "1" and "2" show different bucklings. The Al atoms 
labeled by "A" and "B" are used to define the plane for plots of density profiles as explained in 
the text. 

the hollow site above the substrate, i.e., above the "D" atom shown in Fig. [IJ and has 
occupied the Al-evacuated position. The total energy of the new equilibrium state (state 
"/") was 0.18 eV higher than the state when the Au atom lies on the substrate (state "z"). It 
is therefore, possible that increasing the temperature enhances this and similar mechanisms 
which are responsible for diffusion. We have used the nudged elastic band method 30 to 
determine the minimum energy path for the transition from state i to state /. For this 
purpose, a discretized path consisting of 6 replicas of the system was constructed using 
linear interpolation between the specified i and / states and let to be optimized iteratively. 
In Fig. [7| we have shown the variation of the energy along the minimum energy path. The 
activation energy barrier, E a for this process which is the difference between the maximum 
energy along the path and the energy of the initial state, is obtained to be 0.63 eV. We 
have also used this method to obtain the energy barrier for hopping of Au atom from one 
hollow site at the surface to the neighboring one and obtained a symmetric curve with 
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" ' 0-0 0.2 0.4 0.6 0.8 1.0 

Reaction coordinate 

FIG. 2: Potential energy along the minimum energy path. The left and right endpoints correspond 
to states i and /, respectively. 

E a = 0.64 eV, which is of the same order as that of the latter process. Assuming that 
the diffusion can be described within harmonic transition-state theory,— the transition rate, 
1/r, at a temperature T can be written as 

- = v e ~ E ^ T , (6) 

T 

in which E a is the activation energy barrier, and v is an effective vibrational frequency 
which for diffusion on metal surfaces is of order 10 12 s -1 .— Using the activation energy from 
Fig. [2X we obtain the diffusion rate of ~ 20 s _1 at room temperature, which is negligible; 
and significant value of ~ 10 5 s _1 at 200° C in good agreement with experiment.- 



C. Geometry of adsorbed Au atoms on Al(OOl) surface 

We have considered different supercells with their respective surface unit cells as shown 
in Fig. [H to study the adsorption of the gold atoms on the Al(001) surface. In Table IHH 
we have presented the relaxations of the first and second layers of the Al(001) substrate for 
different coverages. In the case of nonzero coverage, since the atoms in the Al(001) layers 
experience different surroundings, their displacements could be different and therefore, the 
relaxations are calculated for the mean position of the layers. As is seen from Table 1111} for 
0=0.11 ML, the relaxations are still as expansions while, in 0=0.25 ML, the first spacing 
is expanded and the second is contracted. In B=0.5 ML, the first is contracted while the 
second expanded; in 0=0.75 ML, both are contracted; and finally, in 0=1.00 ML, the 
first spacing is contracted, whereas the second one is expanded. To go further from the 
mean behaviors, we mention that the atoms of the first layer of Al(001) lie in a plane for all 
coverages except for 0=0.11 ML, in which they show a small bucklings (<< 0.01 A); and the 



TABLE III: Mean relaxations of the Al(OOl) substrate for different coverages. Here, the average 
of the atomic positions of a layer in the adsorbate-substrate system is taken as the mean position 
of the layer. The = 0.00 ML represents the clean surface. 







9 (ML) 






Aij(%) 0.00 


0.11 


0.25 


0.50 0.75 


1.00 


A 12 +0.91 


+0.68 


+0.25 


-0.21 -0.94 


-1.21 


A 23 +0.26 


+0.15 


-0.43 


+0.59 -0.58 


+0.01 


TABLE IV: Mean distance 


of the Au plane 


from the Al(001) surface for different coveraj 


res. The 


remaining is the same as in 


Table QUI 












G (ML) 






0.00 


0.11 


0.25 


0.50 0.75 


1.00 


d Au -Ai(A) 


1.55 


1.54 


1.67 1.73 


1.80 



second layer shows different behaviors for different coverages. For example, in 0=1.00 ML, 
all the atoms in the second layer lie in a plane; in 0=0.5 ML, the buddings are vanishingly 
small (<< 0.01 A); in 0=0.11 ML, the atoms labeled by "C" lie 0.01 A above the plane and 
the atoms just under the Au atoms lie 0.02 A below the plane while, the displacement of 
the others are negligible. In 0=0.25 ML, the "C" atoms lie 0.02 A above the plane, the "D" 
atoms lie 0.01 A below the plane, and the atoms just under the Au atoms lie 0.03 A below 
the plane. Finally, in 0=0.75 ML, the Al atoms under the "l"-labeled Au atoms, Au(l), 
lie 0.03 A below, the Al atoms under the "2"-labeled Au atoms, Au(2), lie 0.03 A above, 
and the "C" atoms lie 0.04 A above the second Al plane. It is seen from Table HVl that 
the distance of the Au adlayer from the Al substrate almost increases with the coverage. In 
coverage 0=0.75 ML, the Au atoms are grouped in two classes. The Au(l) atoms lie 0.06 A 
below the mean Au adlayer plane while, the Au(2) atoms lie 0.12 A above that plane. 

D. Au adsorption energies for Al(OOl) surface 

Using Eq. ([1]), we have calculated the adsorption energy of gold atoms on the Al(001) 
surface for different surface densities (0.11 ML< <1.00 ML). The result is summarized in 
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TABLE V: Adsorption energies of Au atoms for different coverages. 





0.11 


0.25 


e (ml) 

0.50 


0.75 


1.00 


E b (eV) 


3.72 


3.73 


3.78 


3.77 


3.76 



Table PVT. As is seen, the adsorption energy increases with coverage up to G = 0.5 ML; and 
from 0=0.5 ML to 0=1.00 ML it has a decreasing behavior. However, the differences are 
so tiny (at most 0.06 eV) that practically they do not change with the coverage. 

E. Work functions and surface dipole moments 

Using the self-consistent electrostatic potentials and the Fermi energies, we have calcu- 
lated the work functions, W; and thereby, using the Helmholtz relation [Eq. (jHJ)], we have 
calculated the values of the induced surface dipole moments, //, of the Au/Al(001) surface 
with coverages 0=0.11, 0.25, 0.50, 0.75, 1.00 ML (method 1). The results are summarized 
in Table ED 

As is seen in Table |VT| from 0=0.11 ML to 0=0.25 ML, the change in the work function 
is quite small (+0.02 eV), while, the work-function increases significantly from 0=0.25 ML 
to 0=1.00 ML. On the other hand, the induced dipole moment shows a decreasing behavior 
from =0.11 ML to =0.25 ML with a sharp dip at 0=0.25 ML, an increasing behavior 
from 0=0.25 ML to 0=0.50 ML, and a slow decreasing behavior from 0=0.50 ML to 
0=1.00 ML. Such this "strange" behavior had been also encountered in the earlier work for 
adsorption of oxygen on Cu(lll) by Soon et al. 29 

To proceed in explaining the sharp dip, we have plotted, in Figs. E^a)-(c), the profile of 
the difference electron charge densities in the plane (100) normal to the Al(001) surface, 
which passes through the adsorbed Au, substrate A1(A), and A1(B) atoms. The reason for 
choosing this plane is that the difference electronic charge density was higher than that in 
other planes. As is seen, some features are common in all the three cases. First, there exist 
a region with decrease and a region with increase of charge around the Au atom; second, the 
electronic charges partly reside in the region between the Au atom and the Al substrate; and 
third, there is a decrease of charge around the Al atoms of the first layer. On the other hand, 
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TABLE VI: Work-functions and change in work-functions, in electron-volts; induced surface dipole 
moments, in Debye units, of Au/Al(001) surfaces with different coverages, 0=0.00, 0.11, 0.25, 0.50, 
0.75, 1.00 ML. The 9=0.00 ML refers to the clean surface. 





0.00 


0.11 


V ./ 1 XVXXJ } 

0.25 


0.50 


0.75 


1.00 


w 


4.24 


4.30 


4.32 


4.56 


4.70 


4.81 


AW 


0.00 


0.06 


0.08 


0.32 


0.46 


0.57 


fj, a 




0.13 


0.07 


0.14 


0.13 


0.12 


/"el 




-1.03 


-1.04 


-0.13 


-0.60 


-0.63 


A^ion 




1.15 


1.10 


0.26 


0.72 


0.75 


(Hot 




0.12 


0.06 


0.13 


0.12 


0.12 



a Method 1, calculated using Hclmholtz equation [Eq. ([3|)] 
6 The electronic contribution, the first sum in Eq. ([5]) 
c The ionic contribution, the second sum in Eq. ([5]) 
d The sum of the electronic and ionic parts (method 2) 



the detailed form of the redistributions in different planes normal to the Al(001) plane differ 
from each other. However, since only the ^-component of the dipole moments determines 
the changes in the work-function, we should take an average over the plane parallel to the 
surface. In Figs. IH(a)-(b), we have plotted the averaged difference electron densities as 
functions of z, the normal distance from the central layer of the slab, for different coverage 
structures. To be clear, we have plotted the 6=0.11 ML and 9=0.25 ML cases in Fig. H](a), 
and the remaining cases in Fig. H^b). Inspecting the positions of the adsorbed Au atoms 
(the average position in the case of 9=0.75 ML) in the adsorbate-substrate systems reveals 
that the local maxima to the left of the outermost minima in Figs. H(a)-(b) are located at 
the positions of adsorbed Au atoms. In Fig. H](a), the difference densities at the Au atoms 
are negative (i.e., the Au atom in the adsorbate-substrate system has a less electron density 
compared to the Au atom in the isolated adlayer) while, in Fig. Hl^b), they are positive. 
9n the other hand, in the case of 9=0.25 ML, the value of the difference density at the 
position of the Au atom is minimum compared to the lower and higher coverage structures. 
Using the z-dependent difference electron densities, and multiplying by the slice thickness 



12 




FIG. 3: Contour plots of the difference electron densities in the (100) plane normal to the Al(001) 
surface which passes through the Au atom and its two nearest A1(A) and A1(B) atoms (see Fig.Q]). 
(a), (b), and (c) correspond to 0=0.11, 0.25, and 0.50 ML, respectively. Contour levels in (a) 
and (b) go from -0.03 to 0.05 e/(a.u.) 3 in steps of 0.01 e/(a.u.) 3 ; and in (c), from -0.02 to 0.03 
e/(a.u.) 3 in steps of 0.01 e/(a.u.) 3 , while those between -0.002 and 0.002 e/(a.u.) 3 were taken in 
steps of 0.001 e/(a.u.) 3 . Solid (broken) lines indicate positive (negative) values. The thicker solid 
lines indicate the zero contour levels. The solid circles represent Al atoms. 

(5z = 0.01 A) and the slice position, Zi, we have obtained the electronic contributions, /z e i, 
of the induced surface dipole moments [first term in Eq. (jSJ)]. The electronic and the ionic 
contributions as well as the sum of them are presented in Table I VI I It is clear from Table IVII 
that the electronic contribution has a sharp maximum at G=0.50 ML which can not be 
related to the sharp dip in the induced surface dipole moment. However, adding the ionic 
part yields the correct values of the induced dipole moments, /i to t, which we had called as 
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FIG. 4: Difference electron densities averaged over the surface unit cell; (a)- for 0=0.11 ML and 
0=0.25 ML, and (b)- for 0=0.5, 0.75, 1.00 ML. The positive values imply an increase in the 
electronic density. 

method 2. The results show an excellent agreement with those obtained in method 1 (using 
Helmholtz relation). 



F. Local density of states 

Analyzing the local density of states (LDOS) of Al atoms at the central layer of the 
slab and Al atoms at surface shows some shifts in energies and heights. The situation for 
coverage G = 0.25 ML is shown in Figs. [5] (a)-(b) for s and p orbitals, respectively. The 
Fermi energy of the Au/Al(001) system for coverage 0.25 ML is obtained as 0.08 eV. As 
is seen in Fig. [5] (b), for energies near to Fermi energy, the density of p states shows an 
increase for the Al atoms in the surface layer. On the other hand, in energies much lower 
than the Fermi energy, the p-DOS decrease for the surface Al atoms. However, the s-DOS 
in Fig. [5] (a) does not show much difference at the Fermi surface, but some energy shifts at 
lower energies. Moreover, it is seen that the density of states (both s and p) for the surface 
Al in the Au/Al(001) system are separated in to two groups with a gap around -4.5 eV. 
Conventionally, the lower-energy bunch is called as "bonding" states, and the higher one as 
"anti-bonding" states. 

In order to see the effects of Au adlayer on the surface Al atoms, we have compared, in 
Figs. EKa)-(b), the s- and p-BOS of the Al atoms in the clean Al(001) and the Au/Al(001) 
systems with coverage G=0.25 ML. The plots show that the above-mentioned bunching 
appears because of the adsorbed Au atoms on the surface. The Fermi energy for the clean 
Al(001) system is -0.33 eV, and in the plots we have identified the Fermi energies as the 
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FIG. 5: (a) and (b) represent LDOS for s and p orbitals, respectively, of Al atoms in the bulk 
(central layer of the slab) and surface layer of the Au/Al(001) system with = 0.25 ML. 




Energy (eV) - - . Energy (eV) 

FIG. 6: (a) and (b) represent LDOS for s and p orbitals, respectively, of the surface Al atoms in 
the clean Al(001) and Au/Al(001) surfaces for = 0.25 ML. 

zero of energy. In Fig. E^a), the comparison of the s-DOS's shows that near the Fermi 
energy, there is a significant decrease for the Au/Al(001) system, while for lower energies 
(< —5.5 eV) there is a small increase. On the other hand, in Fig. [6](b) we see that the 
p-DOS of the Au/Al(001) system has a small decrease near the Fermi surface, while it has 
an increase for lower energies (<-4.5 eV). Since, the p-DOS of the Au/Al(001) system has a 
peak at ~ -5 eV, there will be a small increase of the p charge on the Al atom in Au/Al(001) 
system. 

It is also instructive to compare the Au-DOS and the Al-DOS for the nearest neighbor 
Al and Au atoms at surface. We have shown this comparison for the case of G = 0.25 ML 
in Fig. [71(a)- (d). Looking at Figs. [TJa) andUJb) shows a significant overlap between s(Au) 
and p(Al) on the one hand, and a significant overlap between p(Au) and s(Al), on the other 
hand, as shown in Figs. [TJc) and [7J(d). Here we see that the bunching also happens in the 
Au-DOS's. At the same time, the sharp peaks (at ~ -5 eV) of s(Au) and p(Al) and the 
sharp peaks (at ~ -6 eV) of p(Au) and s(Al) indicate a strong hybridization between the 
corresponding two orbitals. 
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FIG. 7: LDOS of the s and p orbitals for surface Al and Au atoms for =0.25 ML. As is seen, 
the sharp peaks of s(Au) and p(Al) indicate a strong hybridization between the two orbitals. The 
Fermi energy is identified as the zero of energy. 

To analyze the bonding of Au in the Au/Al(001) system, we have further compared the 
s-DOS and p-DOS of the Au atom in the isolated Au adlayer with those corresponding to the 
Au atoms in the Au/Al(001) system for 0=0.25 ML. Our results for the isolated Au adlayer 
show that the s-DOS (with large peak) and the p-DOS (with a small peak) are centered 
at 0.16 eV and 0.06 eV with respect to the Fermi energy, and are partially filled. On the 
other hand, Figs. [3a ) and [7(c) show that the electrons completely occupy the bonding levels 
and partially occupy the anti-bonding levels, which lie lower than the energy position of the 
peaks in the isolated adlayer system. This redistribution of the electrons to the lower energy 
levels leads to the bonding of the Au atoms on the Al(001) surface. 

Finally, in Fig. [H] we have depicted the evolution of LDOS for the adsorbed Au atoms 
with different coverages. In this figure, the Fermi energy is identified as the zero of energy. 
For the system with coverage 0=0.75 ML, because of two inequivalent Au atoms "1" and 
"2", we have two curves in the plots which are slightly different in each case. As is seen, 
with increasing O, the density of states undergo shifting and broadening. This broadening 
and shifting is significant for the rf(Au) orbital. 
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FIG. 8: Evolution of the LDOS for the s, p, and d orbitals of Au adsorbed atom with different 
coverages. For coverage 0=0.75 ML, because of two inequivalent Au atoms "1" and "2", we have 
two curves in the plots which are slightly different in each case. As is seen, with increasing 6, the 
density of states undergo shifting and broadening. This broadening and shifting is significant for 
the g?(Au) orbital. 

IV. CONCLUSIONS 

In this work, we have studied the effects of Au adsorption on the surface properties of the 
Al(001), by performing ab initio calculations in the framework of density functional theory 
and supercell methods. Having calculated the adsorption energies, we have found the hollow 
sites on the surface as the most preferred adsorption sites, and have considered the adsorption 
of Au atoms at the hollow sites for different coverages from = 0.11 ML to = 1.00 ML. 
We have shown that the adsorbed Au atoms, because of higher electronegativity, attracts 
the electronic charge of the substrate and cause the induction of surface dipole moment. 
The induction of the surface dipole moment gives rise to the work-function change. Using 
the Helmholtz relation, we have obtained a sharp dip in the induced dipole moment at 
0=0.25 ML. We have shown that using the self-consistent electronic density and taking into 
account the ionic charge redistributions, one can directly calculate the induced surface dipole 
moment. Therefore, the sharp dip is not only due to electronic charge redistributions but, 
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one should take into account the ionic redistributions as well. Moreover, we have shown that, 
at zero temperature, the Au atoms do not diffuse into the Al(OOl) substrate, in agreement 
with the experimental observations. To gain an insight on the nature of bondings, we have 
also studied the local density of states. Finally, for sufficiently high Au coverages, one has 
to consider the clustering effects. The clustering leads to the formation of Au nano-clusters 
on the Al(OOl) surface which is important in catalytic applications. Work in this direction 
is in progress. 
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